ASYMPTOTICS FOR DISCRETE WEIGHTED MINIMAL 
RIESZ ENERGY PROBLEMS ON RECTIFIABLE SETS 



S.V. BORODACHOV\ D.P. HARDIN^ AND E. B. SAFF^ 

Abstract. Given a compact d-rectifiable set A embedded in Euclidean 
space and a distribution p{x) with respect to d-dimensional Hausdorff 
measure on A, we address the foUowing question: how can one gener- 
ate optimal configurations of A'^ points on A that are "well-separated" 
and have asymptotic distribution p{x) as N ~> oo? For this purpose 
we investigate minimal weighted Riesz energy points, that is, points in- 
teracting via the weighted power law potential V — ■w{x,y) \x — J/l"", 
where s > is a fixed parameter and w is suitably chosen. In the un- 
weighted case {w = 1) such points for A'^ fixed tend to the solution of 
the best-packing problem on A as the parameter s oo. 



1. Introduction. 

Points on a compact set A that minimize certain energy functions often 
have desirable properties that reflect special features of A. For A = S'^, the 
unit sphere in M^, the determination of minimum Coulomb energy points is 
the classic problem of Thomson j2L)l [3] . Other energy functions on higher 
dimensional spheres give rise to equilibrium points that are useful for a vari- 
ety of applications including coding theory jH], cubature formulas f21_, and 
the generation of finite normalized tight frames In this paper, we shall 
consider a generalized Thomson problem, namely minimum energy points 
for weighted Riesz potentials on rectifiable sets. Our focus is on the hyper- 
singular case when short range interactions between points is the dominant 
effect. Such energy functions are not treatable with classical potential the- 
oretic methods, and so require different techniques of analysis. 

Let ^ be a compact set in M*^ whose d-dimensional Hausdorff mea- 
sure, 7id{A), is finite. For a collection of A^(> 2) distinct points ujm := 

1991 Mathematics Subject Classification. Primary 11K41, 70F10, 28A78; Secondary 
78A30, 52A40. 

Key words and phrases. Minimal discrete Riesz energy. Best-packing, Hausdorff mea- 
sure, Rectifiable sets, non-Uniform distribution of points. Power law potential. Separation 
radius. 

^Research of this author was conducted as a graduate student under the supervision of 
E.B. Saff and D. P. Hardin at Vanderbilt University. 

^The research of this author was supported, in part, by the U. S. National Science Foun- 
dation under grants DMS-0505756 and DMS-0532154. 

^The research of this author was supported, in part, by the U. S. National Science Foun- 
dation under grant DMS-0532154. 

1 



2 



S.V. BORODACHOV, D.P. HARDIN, AND E. B. SAFF 



{xi, . . . jXn} G a, a non-negative weight function w on A x A (we shall 
specify additional conditions on w shortly), and s > 0, the weighted Riesz 
s-energy of ujn is defined by 



where \X\ denotes the cardinality of a set X. Since, for the weight w{x, y) := 
{w{x,y) + w{y, x))/2, we have 



we shall assume, without loss of generality, throughout this paper that w is 
symmetric, i.e., w{x,y) = w{y,x) for x,y £ A. We call w : A x A ^ [0,oo] 
a CPD-weight function on ^ x ^ if 

(a) w is continuous (as a function on yl x ^) at "H^-almost every point 
of the diagonal D{A) := {{x,x) : x G A}, 

(b) there is some neighborhood G of D{A) (relative to Ax A) such that 
infc w > 0, and 

(c) w is bounded on any closed subset BcAxA such that Br\D{A) = 0. 

Here CPD stands for (almost) continuous and positive on the diagonal. 
In particular, conditions (a), (b), and (c) hold if w is bounded on A x A 
and continuous and positive at every point of the diagonal D{A) (where 
continuity at a diagonal point {xq, xq) is meant in the sense of limits taken 
on ^ X A). 

If w = 1 on A X A (which we refer to as the unweighted case), we write 
Es{ujn) and £s{A, N) for Ef{u;N) and £f{A, N), respectively. For the trivial 
cases = or 1 we put Es{ujn) = SsiA,N) = Ef{uJN) = £^{A,N) = 0. 

We are interested in the geometrical properties of optimal s-energy 
A^-point configurations for a set A; that is, sets wtv for which the infimum 
in is attained. Indeed, these configurations are useful in statistical 
sampling, weighted quadrature, and computer-aided geometric design 
where the selection of a "good" finite (but possibly large) collection of 
points is required to represent a set or manifold A. Since the exact 
determination of optimal configurations seems, except in a handful of cases 
(cf. [221 El Q El 121 ini) 1 beyond the realm of possibility, our focus is on the 
asymptotics of such configurations. Specifically, we consider the following 
questions. 

(i) What is the asymptotic behavior of the quantity £^{A, N) as N gets 
large? 

(ii) How are optimal point configurations distributed as ^ oo? 




while the N -point weighted Riesz s-energy of A is defined by 
(1) SfiA, N) := mi{Ef{iON) : ujn C A, \ujn\ = N}, 



Ef{u;N) = Ef{u;N) = 2 Yl 



l<i<j<N 



w{Xi,Xj) 

I I ^ 

X i X ^ 
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(iii) What estimates can be given for the minimal pairwise distance 
between points in optimal configurations for large A^? 

In the unweighted case, much is known regarding these questions. In 
particular, when s < dim A (the Hausdorff dimension of A), the limit distri- 
bution of optimal A^-point configurations is given by the equilibrium measure 
Xs^A that minimizes the continuous energy integral 

hil^) ■=[[ I — diJ.{x)dn{y) 
J J \x-y\ 

AxA 

over the class M.{A) of (Radon) probability measures fi supported on A. In 
addition, the asymptotic order of the Riesz s-energy is N"^; more precisely 
we have £s{A,N)/N^ Is{Xs,a) as iV ^ oo (cf. |IHI Section II.3.12]). In 
the case when A = 5"^, the unit sphere in M'^"'"^, the equilibrium measure is 
simply the normalized surface area measure. 

If s > dim A, then IsifJ-) = oo for every /i S A4{A) and potential theoretic 
methods cannot be used. However, by using techniques from geometric 
measure theory, it was recently shown in J2] that when A is a d-rectifiable 
manifold of positive d-dimensional Hausdorff measure and s > d, optimal 
A^-point configurations are uniformly distributed (as oo) on A with 

respect to d-dimensional Hausdorff measure restricted to A. The assertion 
for the case s = d further requires that ^4 be a subset of a manifold (see 
Theorem E] in Section 1.1). 

Our motivation for considering the weighted minimal Riesz energy prob- 
lem is for the purpose of obtaining point sets that are distributed accord- 
ing to a specified non-uniform density such as might be used as nodes 
for weighted integration or in computer modeling of surfaces where more 
points are needed in regions with higher curvature. In this paper we shall 
show that for a compact d-rectifiable set A having positive d-dimensional 
Hausdorff measure, A^-point configurations for A minimizing the weighted 
Riesz s-energy are distributed asymptotically with density proportional to 
{w{x, x))~'^/^ provided s > d. (This continues to be true even when w has a 
finite number of zeros on the diagonal, provided their order is less than s.) 
Even in the unweighted case, these results extend those of obtained for 
the class of d-rectifiable manifolds to the more general class of d-rectifiable 
sets. 

For the remainder of this introduction we provide the needed notation 
and discuss known results. Section 2 is devoted to the statements of the 
main results of this paper. The detailed proofs of these main results, which 
utilize the basic lemmas described in Section 3, are provided in Sections 4, 
5 and 6. 

1.1. Notation and previous results. It is helpful to keep in mind that 
minimal discrete s-energy problems can be considered as a bridge between 
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logarithmic energy problems and best-packing ones. Indeed, in the un- 
weighted case {w{x, y) = 1) when s — > and N is fixed, the minimal energy 
problem turns into the problem for the logarithmic potential energy 

Y] log, 1, 

l<i^j<N ' 

which is minimized over all iV-point configurations {xi, . . . ,xm} C A. This 
problem is equivalent to the maximization of the product 



n 



I X 2 j I * 

For planar sets, such optimal points are known as Fekete points. (For the 
case when A = S'^, the polynomial time generation of "nearly optimal" 
points for the logarithmic energy is the focus of one of S. Smale's "problems 
for the next century"; see p^.) 

On the other hand, when s — > oo, and is fixed in the unweighted case, 
we get the best-packing problem (cf. [Hj); i.e., the problem of finding 

A^-point configurations ujn C A with the largest separation radius: 

(2) 6{ujn) '■= min \xi — Xj\. 

In this paper we will consider the case s > dim A. Let C^' be the Lebesgue 
measure in M*^ and Tid be the d-dimensional Hausdorff measure in M'^ nor- 
malized so that its restriction to M*^ C M*^ is Cd- Denote by Bd'{xo,r) the 
open ball in M"^ centered at the point xq with radius r > and set 

(3) ^^-^^™0'^)) = ZIW)- 

Given sequences {a7v}^=i and {6Ar}^^]^ of positive numbers, we will write 
ajv ~ b]\[, N oo, if limTv^oo ciN/bN = 1- 

Regarding questions (i) and (ii), A.B.J. Kuijlaars and E.B. Saff |E] proved 
that for the unit sphere S'^, there holds 

£diS^N)r.-J±^^N'logN, N^oo, 

and it is known that the distribution of the minimal energy points is asymp- 
totically uniform in this case. For one-dimensional rectifiable curves in R'^ , 
the paper by A. Martinez-Finkelshtein et al [T^ provides answers in the 
unweighted case to questions (i) and (ii), as well as question (iii) for regular 
Jordan arcs or curves. 

Question (iii) in the unweighted case has also been considered for several 
other special cases. B.E.J. Dahlberg [7] proved that if ^ = 5"^ and s = d—1, 
d > 2, or ^ C is a smooth surface and s = 1, d = 2, there is a constant 
C > such that for every s-optimal collection C A with A^ points 

(4) 6{uj*n) > CN-^/^. 
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In it was shown that (HJ holds for A = S'^ when s > d. 

In J2] and ^Hl, questions (i), (ii), and (iii) were addressed for a more 
general class of sets A which we now describe. First recall that a mapping 
(/) : T — > R'^ , T C K'^, is said to be a Lipschitz mapping on T if there is some 
constant A such that 

(5) |(/)(x) - (l){y)\ < \\x - y\ for x,y eT, 

and that (p is said to be a bi-Lipschitz mapping on T (with constant if 

(6) {l/X)\x-y\<\(t>ix)-(t>iy)\<Mx-y\ forx,yGT. 

Following |12_, we say that a set A C M^' is a d-rectifiable manifold if it is a 
compact subset of a finite union of bi-Lipschitz images of open sets in M'^. 

We further recall that if A C M'^ is compact and u and {i^n}'n=i Borel 
probability measures on A, then the sequence converges weak-star to ly 
(and we write z/jv * > i^) if for any function / continuous on A, we have 

lim / fdi^N = / fdiy. 

Denote by 6x the atomic probability measure in M'^' centered at the point 
X e R'^'. 

For future reference and comparison we now state some main results from 

m- 

Theorem A. Let A C M^' be a compact d-rectifiable manifold and suppose 
s > d. Then 

where Cs^d is a positive constant independent of A. 

Furthermore, if 7id{A) > 0, any asymptotically s-energy minimizing se- 
quence of configurations ujn = {^i j ■ ■ ■ i x^}, ^ = 2, 3, . . ., for A is uni- 
formly distributed with respect to TLd; that is, 

1 ^ 



'Hd\A 

N > , . , , iV oo. 



By asymptotically s-energy minimizing we mean 

Remark. For s > d, the constant Cg^d appearing in Q of Theorem ^ 
can be represented using the energy for the unit cube in via formula Q: 

For d = 1 and s > 1, it was shown in ^\ that Cs,i = 2C(s), where C(s) 
is the classical Riemann zeta function. However, for other values of d, the 
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constant Cs^ is as yet unknown. For the case d = 2, it is a consequence of 
results in ITO that 



(9) Cs,2<{V^/2X'\l{s), 



where Cl{s) is the zeta function for the planar triangular lattice L consisting 
of points of the form m(l, 0) + n(l/2, \/3/2) for m,n ^ Z. It is conjectured 
in US] that in fact equality holds in ©. Furthermore, it is shown in Uj that 

as s — > oo 

iCs,2f' ^ V3/2, 

which is consistent with this conjecture. 

When < 7id{A) < oo we observe that the minimum energy experiences 
a transition in order of growth; namely, as s increases from values less than d 
to values greater than d the energy switches from order A^^ to order A^^+'^Z'^ 
as — > oo. As the following theorem from 12 describes, at the transition 
value s = d, the order of growth is A^^logA^. In the proof of this fact, a 
more delicate analysis was utilized that involved an additional regularity 
assumption on A. 



Theorem B. Let A be a compact subset of a d- dimensional -manifold in 
W^' . Then 

(10) lim = -li-, 

where [3d is the volume of the d- dimensional unit ball as defined in 

Furthermore, if TCd{A) > 0, any asymptotically d-energy minimizing se- 
quence of configurations ujn = {x^ , ■ ■ ■ , x^}, N = 2, 3, . . for A is uni- 
formly distributed with respect to 7id; that is, ^ holds. 

Regarding separation results, the following was shown in 12 . If A C 
M.'^ is a bi-Lipschitz image of a compact set from M'^ of positive Lebesgue 
measure, then for every s > d there is a constant > such that 



'csN-^/'^, s>d, 
CdiN\ogN)-^/'^, s = d, 



for every s-optimal A^-point configuration on A. This result was extended 
by S. Damelin and V. Maymeskul in jSj to a finite union of bi-Lipschitz 
images of compact sets from M'^ . 



2. Main results. 

In this paper we extend Theorem A to the class of d-rectifiable sets; 
where, by a d-rectifiable set ^ C M'^ , we mean the image of a bounded set 
in under a Lipschitz mapping (cf. JHI)- Consequently, we relax the 
bi-Lipschitz condition in Theorem A. Furthermore, both Theorems A and 
B are extended to the case of weighted energy, and separation estimates for 
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optimal configurations are obtained even when the Hausdorff dimension of 
the compact set is not necessarily integer. 

Our first main result extends of Theorem A to the class of d-rectifiable 
sets. 

Theorem 1. Suppose A C M'^' is a closed d-rectifiable set. Then, for s > d, 

where Cg^d is as in Theorem A. 

A key component of this extension is a regularity result given in 
LemmalSJof Section 3 showing that liminfAr^oo£^s(i(', A^)/A/'^"'"'*/'^ is close to 
liminfAf^oo-^sC^, A^)/A^-^+''/'^ for compact sets K C A such that HdiA \ K) 
is sufficiently small. 

If ^ is a compact set in and w is a CPD-weight function on A x ^, 
then for s > dwe define the weighted Hausdorff measure TC''^^ on Borel sets 
BcAhy 

(12) n'^^{B):= [ {w{x,x))-''/'dnd{x), 

Jb 

and its normalized form 

(13) h'^'"iB) ■.= n'j'"iB)/n'^'"{A). 

We say, that a sequence {ujN}'N=i of A^-point configurations in A is asymp- 
totically (w, s)-energy minimizing for A if 

E^iuN)^£:iA,N), N^oo. 

The main results of this paper include the following generalizations of 
Theorems EI and El 

Theorem 2. Let A C R'^' be a compact d-rectifiable set. Suppose s > d and 
that w is a CPD-weight function on Ax A. Then 

(14) lim 



where Cs,d is as in Theorem rM 

Furthermore, if Ti.d{A) > 0, any asymptotically {w, s)-energy minimiz- 
ing sequence of configurations ujn = {x^ , . . . , x^}, N = 2, 3, . . ., for A is 
uniformly distributed with respect to 'H^/" ' ^^^^ 



N 

* 



(15) 1^^^,^/,^, AT^oo. 



^^-x^- "d ' 

k=l 
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Theorem 3. Let A be a compact subset of a d- dimensional -manifold in 
and suppose w is a CPD-weight function on A x A. Then 

(16) lim ^^"('^•"> - 



N^oo N'^logN n'^''"{A)' 

Furthermore, if TC^{A) > 0, any asymptotically {w,d)-energy minimiz- 
ing sequence of configurations ujn = {x^ , . . . , x^}, N = 2,3, . . ., for A is 
uniformly distributed with respect to Ti.'^^ ; that is, hl!^) holds with s = d. 

Remarks. In the case 7i(i{A) = 0, the right-hand sides of H14|) and (|16() 
are understood to be infinity. 

Next we obtain estimates for the separation radius of optimal configu- 
rations on sets of arbitrary Hausdorff dimension a. We remark that the 
normalization for the Hausdorff measure 7ia plays no essential role here. 

Theorem 4. Let < a < d' . Suppose A C R'^ is a compact set with 
7ia{A) > and let w be a CPD-weight function that is bounded and lower 
semi- continuous on A x A. Then, for every s > a there is a constant 
Cs = Cs{A,w,a) > such that any {w, s)-energy minimizing configuration 
:= {xi^N, ■ ■ ■ ,xn,n} on A satisfies the inequality 

* , - I I ^ jcsN-y^, s > a, 

OlWiv = mm \Xi N — Xj m\ > < , , , , 

i<i^j<N^ ' ' - |c«(A^logiV)-V", s = a, iV>2. 

As a consequence of the proof of Theorem 0] we establish the following 
estimates. Let 

n^{A) := inf{^ (diam G^)" : ^ C (J d}. 

i i 

Corollary 1. Under the assumptions of Theorem^ for N > 2, 

' Ms,a\\w\\AxAn'^{A)-'/''N^+-'/", s > a. 



M„N^los.N, s = a. 



where the constant Mg^a > is independent of A, w and N , and the constant 
Ma is independent of N . 

In order to obtain a finite collection of points distributed with a given 
density p{x) on a d-rectifiable set A, we can take any s > d and the weight 

(17) wix, y) := (p(x)p(y) + \x - y\r'/^'', 

where the term |x — y| is included to ensure that w is locally bounded 
off of D{A). By Theorems |2l and 01 any asymptotically (it;, s)-energy 
minimizing sequence of A'^-point configurations will converge to the required 
distribution as N ^ oo. We thus obtain 
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Corollary 2. Let A C M^' be a compact d-rectifiable set with TCd{A) > 0. 
Suppose p is a bounded probability density on A (with respect to TLd) that is 
continuous TC^- almost everywhere on A. Then, for s > d and w given by 
jlT^ , the normalized counting measures for any asymptotically (w, s)-energy 
minimizing sequence of configurations ojm converge weak* (as N ^ oo) to 
pdTCd- 

Furthermore, if inf^ p > and p is upper semi- continuous, then any 
{w, s)-energy minimizing sequence of configurations ojn is well- separated in 
the sense of Theorem 4 with a = d. 

Remark: The first part of Corollary|2lholds for s = d when A is contained 
in a d-dimensional manifold. 

Finally, we consider weight functions with isolated zeros. For t > 0, we 
say that a function w : A x A ^ M. has a zero at (a, a) G DiA) of order at 
most t if there are positive constants C and 5 such that 

(18) w{x,y) > C\x - a\^ {x,y e An Bd'{a,5)). 

If w has a zero a £ A whose order is too large, then a may act as an attractive 
"sink" with £'^{A,N) = 0. For example, let A be the closed unit ball in 
W^, vu{x, y) = + for x,y £ A with t > s > d. If lon = {xi, ■ ■ ■ ,xn) 
is a configuration of points in A, then Ef{'yuj]\f) = 7*~*ii^^(a;Ar) for any 
< 7 < 1. Taking 7^0, shows that £^{A, N) = 0. 

A closed set A C M'^ is a-regular at a £ A if there are positive constants 
Co and 5 such that 

(19) (Co)-V" < nM^Bd'{x,r)) < C^r^ 
for all x G ^ n B^; (a, b) and < r < (5. 

Theorem 5. Let A C M'^' be a compact d-rectifiable set and s > d. Suppose 
A is ai-regular with ai < d at ai, i = 1, . . . ,n, for a finite collection of points 
ai, . . . ,an in A and that w : A x A ^ [0,oo] is a CPD-weight function on 
K X K for any compact K d A \ {ai, . . . , a„}. // w has a zero of order at 
most t < s at each {ai,ai), then the conclusions of Theorem\^ hold. 

Remark: The hypotheses of Theorem El imply that 
{w{x,x))^^/^dHd{x) < oo (see Section 6). 

3. Lemmas. 

In this section we prove several lemmas which are central to the proofs of 
our main theorems. 

3.1. Divide and conquer. In this subsection we provide two lemmas re- 
lating the minimal energy problem on A = B U D to the minimal energy 
problems on B and D, respectively. 
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In order to unify our computations for the cases s > d and s = d, we 
define, for integers > 1, 



N^+'^l'^, s> d, 
log N, s = d 



and set Tg^iN) = 1 for = or 1. For a set A C and s > d, let 
and 

if this hmit exists (these quantities are ahowed to be infinite). In the case 
w{x,y) = 1, we use the notations ^(^)) 5s,d(^) ^^'^ gs,d{^)i respectively. 

Let dist(i?,Z?) := inf {|x — y\ : x £ B, y £ D} denote the distance be- 
tween sets B,D C M"^ . The following two lemmas extend Lemmas 3.2 
and 3.3 from |12j to the weighted case. We remark that the following re- 
sults hold when quantities are or infinite using 

Lemma 1. Let s > d > and suppose that B and D are sets in M'*' such 
that dist(i?, D) > 0. Suppose w : {B U D) x {B U D) ^ [0, oo] is bounded on 
the subset B x D. Then 

(20) gl,{B U D)-'/^ > glABY"'' + gldiDY"''. 

Proof. Assume that < g'^^B), g'^^iD) < oo. Denote 



a 



g^/BY/^+g^^DY/s- 



For G N, let Nb '■= [a*N\ (where [x\ denotes the greatest integer less 
than or equal to x), No := N — Nb and C B and uj^ C D he configura- 
tions of Nb and Nd points respectively such that Ef{u}^) < £f{B, Nb) + 1 
and Ef{LO^) < £f{D,ND) + 1. Let 70 := dist(S,L>) > 0. Then 

£1^{B^D,N)<E^{lo^Uu:^) 

^-^ \x — y\ 

< £fiB,NB)+£T{D,ND)+2 + 2%'N^\\w\\BxD, 

where ||bxD denotes the supremum of w over B x D. Dividing by Tg^N) 
and taking into account that TsANb)/tsAN) — > {a*Y~^'^^'^ as N —>■ 00, we 
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obtain 



5s,d(^UL>) <limsup — — +limsup — — 

= lim sup ■ + lim sup — - — 

—i-oo '^s,d{-^^ ) 

< gUB) . {a*)'"-'" + 91AD) ■ (1 - a*)'-^-^/" 



-d/s 



-s/d 



The remaining cases when g^^B) or g^^D) are or oo easily follow from 
the monotonicity of g^^- D 

The following statement in particular shows sub-additivity of 9^ ■ 

Lemma 2. Let s > d > and B,D C R'^' . Suppose w : (BUD) x (BUD) 
[0,oo]. Then 

(21) gl^B U D)-'/^ < fjB)-'^'' + gl.iD)-"'. 

Furthermore, if g"^ ^B) ^ g"^ ^D) > and at least one of these quantities is 
finite, then 

(22 lim ^— ^ = , 

Af^N^oo N g^JyBYl'' + g-^JyDYi" 

holds for any sequence {uiN^NeN of N -point configurations in B L) D such 
that 

(23) lim ^^=(9:AB)-'^'+9:AD)-'^'^''^' 



where N is some infinite subset ofN. 

In the case g^^D) = oo the right-hand side of relation is understood 
to be 1. 

Proof. Assume that < g"^ J^B) , g"^ J^D) < oo (we leave other cases to the 
reader). Let an infinite subset A/i C N and a sequence of point configurations 
{LON}N€Afn C B U D, he such that lim^^gTv^oo I'^A'' H B\/N = a, where 
< a < 1. Set Nb := \ujn n B\ and No ■= \u;n \ B\. Then 

EfM > EfiuJN nB) + EfiuJN \ B) > £f{B, Nb) + £T{D, Nd), 

and we have 



MlBN^cx. TsA^) M97V^«D Ts^Nb) Ts^N) 
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(24) > F{a) := g;^^^{B)a'-^^/Ug:^/D){l-ay+^/''. 

Let 



a 



and {uiN^NeU be any sequence of point sets satisfying (pS)) . If J\f2 C AA is 
any infinite subsequence such that the quantity IlD^v r\ B \ /N has a hmit as 
J\f2 3 N ^ OO (denote it by ai), then in the case ^{B) , ^{D) < oo by 

(ESI) and (EU we have 

F(a)= hm ^l%l>F{a^). 

It is not difficult to see that a is the only minimum point of F{t) on [0, 1]. 
Hence ai = 5, which proves 1)22^ . 

Now let {T^N^NeU-i be a sequence of A^-point configurations \n B\JD such 
that 

g-{BiJD)= hm 

(ujn^s can be chosen for example so that Ef{tJiy) < £f{B U D,N) + 1). If 
Mi C A/3 is such an infinite set that lim_v'49Ar^oo l^iv H B\ /N exists (denote 
it by 02), then by (Pl)) we obtain 

> F{a) = {fjB)-'''^+gllD)-^I^Y''\ 

which implies H21() . □ 

3.2. Lemmas from geometric measure theory. Recall that denotes 
the volume of the unit ball in W^. For convenience, we also define /3o := 1- 
For a set C W^' and /i > 0, we let 

W{h) := {x G M"'' : dist(x, W) < h}. 

The upper and the lower Minkowski contents of the set W are defined, 
respectively, by 

r^0+ Pd'-d I ■r^0+ Pd'-d ' 

If the upper and the lower Minkowski contents of the set W coincide, then 
this common value, denoted by A/1^(VK), is called the Minkowski content of 
W. We shall rely on the following property of closed d-rectifiable sets. 

Lemma 3. (see [101 Theorem 3.2.39] IfWc M'*' is a closed d-rectifiable 
set, thenMdiW) = ndiW). 



g'^^iB U Z?) = _ ^ ]im J \ J/ > ^(02) 
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We shall also need the following fundamental lemma from geometric mea- 
sure theory. 

Lemma 4. (see ^TDl, Lemma 3.2.18]j. Let W C M"^' be a d-rectifiable set. 
Then for every e > there exist compact sets Ki, K2, K3, . . . C M*^ and 
bi-Lipschitz mappings ipi : Ki ^ with constant 1 + e, i = 1, 2, 3, . . such 
that ipi{Ki), ip2{K2), -03 (-^^3) J • • • are disjoint subsets of W with 

-Hd (w\\jMKi)^ =0. 

In fact, the above lemma holds for any set of finite W^-measure that is, 
up to a set of "H^-measure zero, the countable union of d-rectifiable sets. 
However, Lemma IHl does not hold for this larger class (cf. 10 ). 



3.3. Regularity Lemma. To get an estimate from below for ^(^) 
will need the following result. 

Lemma 5. Let s > d and suppose A C M'^' is a compact set such that 
Md{A) exists and is finite. Then for every e £ (0, 1) there is some S > 
such that for any compact set K C A with Ad^jK) > j^.d{A) — 5 we have 

(25) g_^/A)>{l-e)g_^^^{K). 

Proof. The assertion of the lemma holds trivially if g^ ^(^) = 00. Hence, we 
assume g^ ^(^) < 00. Let C N be an infinite subset such that 

lim — — 73- = q ,iA). 

Choose 6 £ (0, 1/2^°') and set 

(26) p:=5i/(4d) hN ■.= -p^N-^l\ NeAf. 

3 

Suppose K is a compact subset of A such that M.^(K) > M.d{A) — 6. Then 
there is some Ns G N such that for any N > Ng, N £ J\f, we have 

(27) -j^<Md{A)+5 and > Md{A) - 5. 

For N G M with N > Ng, let cj^ := {xi^n, ■ ■ ■ , xn,n} be an s-energy 
minimizing A^-point configuration on A. For i = 1,...,N, let r-^ : = 
imii\x j^N — Xi^N\ denote the distance from Xi^N to its nearest neighbor in 

io*j^. Further, we partition a;^ into a "well-separated" subset 

and its complement ujj^ := a;^ \ uj^. We next show that uij^ has sufficiently 
many points. (In the case 7id{A) > 0, Theorem |3] implies that |wj\r| = N for 
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sufficiently small p. However, the case 7id{A) = still requires consideration 
and we provide a proof that works in either case). For N G N , we obtain 

N N ^1 

1=1 J=l ' ' Jl I I— I \ I / 

Let ko := ^{A) + 1. There is TVi e N such that for any N > Ni, N e Af, 

For the rest of the proof of this lemma, let G be greater than := 
maxjA'^i, A'^}. Then, 

l^^l ^ 8s{A,N) 
psN - Ni+s/d ^ '^O' 

and, hence, we have 

(28) |wjv| < kop'^N and |cjjv| > (1 - kop'')N. 
Next we consider 

uj% := ujjff^K{3hN), 0J% := u;lf\ K{3hN), 
and show that the cardinality of coj^ is sufficiently large. From ()27() we get 

(29) Cd' [A{hN)\K{hN)] = Cd' [A{hN)] - Cd' [K{hN)] 

< {MdiA) + S)^d'^dhi~'' - (MdiA) - 6)Pd'-dhf-'' 



2(3d'-dSh%-''. 



Note, that 

(30) Fn:= U Bd'ix,hN)cAihN)\KihN). 



',2 



For any distinct points Xi^^^Xj^M £ have 

\xLN - Xj^N\ > rf > pN-^/'^ > p^AT-i/rf = 

Hence, Bdi{xi^iy,hiy)f]Bd'{xj^N,hN) = 0- Then, using (|29|) and (|3Ujl . we 
get 

P%\ = {pd'hiy^ ^d'[Bd'{x,hN)] = [pd'hiy^d'iFN) 

< [pd'hf^y^ Cd'[A{hN)\K{hN)] < 2Pd'-dPd'^5h-N- 
Hence, recalling ()26() . we have 

(31) pl\<2-3''Pd'-dfid'^5^'^N. 
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Let xo := 2 • S'^Pd'-dP^'^ ■ Then, using ^ and (EU, we have 
\^%\ = \^h\- P%\ > (l - kop' - xo5'/') N. 

Next, we choose a configuration of points in K which is close to ujf^ 
and has the same number of points and order of the minimal s-energy as w^. 
For every Xi^jy G pick a point Ui^jy G K such that \xi^]y — yi,N\ < Shjy = 
p2jY-i/rf g^jj(j ._ i^y. j^ ■ g w^}. Since every point Xj^Ar E l<j^ lies 

in wjy, we have 

\xi,N - ViM < p^N~^^'^ < prf < p \xi^N - Xj^nI , j / «• 
Then, if Xi^jy 7^ Xj^jy are points from w^, we have 

\yi,N - yj,N\ = \yi,N - Xi^N + Xi^N - Xj^n + Xj^N - VjM 

> \Xi^N - Xj^nI - \Xi,N - yi,N\ - \xj,N - yj,N\ 

> \Xi^N - Xj^nI - '^p\Xi,N - Xj^nI = (1 - 2p) \Xi^N " Xj^nI ■ 

Since p G (0,1/2), it follows that = \^^%\ and 



T" — 111 



X -y\ ^ \x-y 



>{l-2py Yl ^^ = {l-2prE,{u:%). 



\x - y\ 



Now suppose e S (0,1). We may choose 5 > sufficiently small (recall 
p = 5i/(4d)) so that (1 - 2p)^(l - kop' - xo'5^/^)^+'/'^ > (1 - e)- Hence, 



>(l-2p)^ hminf ^^(^'I'^^l^ 




M3N^oo \, ,2 \'^+s/d 

\^n\ 
> (1 - ^)lJK) 

holds for any compact subset K C A such that M^{K) > Aid{A) — 5. □ 

4. Proofs of Main Theorems. 

Proof of Theorem^ First we remark that ii K CR'^ is compact, then K is 
trivially a d-rectifiable manifold (or set) and so Theorem A shows, for s > d, 



(32) 
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Suppose < e < 1. Since A C M'^' is a compact d-rectifiable set, Lemma^ 
implies the existence of compact sets Ki, K2, i^s, . . . C M'^ and bi-Lipschitz 
mappings : — > R"' , i = 1, 2, 3, . . ., with constant 1+e such that t(^i{Ki), 
'4'2[K2)i '^'^{Ks), . . . are disjoint subsets of A whose union covers "H^-almost 
ah of A. 

Let n be large enough so that 

(n \ n 

\JMK^)] =Y,'>^d{MK^)) > {l + e)-''nd{A). 
i=l J i=l 

Since ipi is bi-Lipschitz with constant (1 + e) we have 

(33) gs,MiKi)) < (1 + eYgsAKi) = CsA^ + eyCdiKi)-'/" 

Applying Lemma Q we obtain 

/ n \ / " \ ^^^'^ 

(34) g.^A) < U ^^(^^) ^ T.dsA^Ki))-''/^ 



\i=l / \i=l 



~s/d 



<C,,rf(l + e)2^ \^^ndiMKi))j <Cs,d'^ + ef'nd{Ay 

We next provide a lower bound for ^(^)- Since A is a closed d-rectifiable 
set, we have M.d{A) = Tlii{A) < 00 (cf. Lemma 01). Let (5 > be as in 
Lemma 121 i.e., inequality H25() holds for every compact set K C A such that 
M-d{K) > M-d{A) — S. Now let n' be large enough so that 

n' 

Md (ytiMKi)) =Y,'^^d[MKi)] > HdiA) - 5 = Md{A) - 5. 

i=l 

As in (|33|) we have 

(35) l,/i^i{Ki)) > (1 + e)-^5s,d(i^i) = C.^l + er^CdiKi)-''" 

yCs^dil + er^'TidrnKi))'''"- 
Then Lemma |^ and H35|) give 
(36) 

Letting e go to zero in ()34() and ()36|). we obtain □ 
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4.1. Proofs of Theorems [21 and ^ The following lemma relates the 
weighted minimal energy problem (s > d) on a set A C M'^ to the un- 
weighted minimal energy problem on compact subsets of A. Theorems El 
and O then follow easily from this lemma. 
For convenience, we denote 

Cd,d-=Pd, den. 

Lemma 6. Suppose s > d, A C M"^' is compact with TCd{A) < oo, and that 
w is a CPD-weight function on A x A. Furthermore, suppose that for any 
compact subset K C A, the limit gs^d{K) exists and is given by 

(37) 9.,diK) = 
Then 

(a) g'^^{A) exists and is given by 

(38) gf^,{A)=Cs,d{nriA))-'/\ 
and, 

(b) if a sequence {ujN}'N=2' '"^^ere ujn = {x^ , • • • > ^n}' ^-^ asymptotically 
(w, s)-energy minimizing on the set A and 7ici{A) > 0, then 

N 



k=l 

Remark. If 7id{K) = 0, condition (|S7)) is understood as gs,diK) = oo. 

Proof. To prove the first part of the theorem, we break A into disjoint 
"pieces" of small diameter and estimate the (t/;, s)-energy of A by replac- 
ing w with its supremum or infinum on each of the "pieces" and applying 
Lemmas ^ and (21 

For 5 > 0, suppose that Vs is a partition of A such that diamP < 5 and 
Hd{P) = TCd{P) for P G Vs, where B denotes the closure of a set B. For 
each P Vs, choose a closed subset Qp C P so that Qs := {Qp : P G Vs} 
satisfies 

(40) '^diQp) > -HdiA) - 

PePs 

An example of such systems Vs and Qs can be constructed as follows. Let 
Gj [t] be the hyperplane in consisting of all points whose j-th coordinate 
equals t. If (— a, a)*^' is a cube containing A, then for i = {ii, . . . € 
{1, . . .,mY', let 

where m and partitions —a = t^ < t\ < ... < tj^ = a, j = 1, . . . , d' , are 
chosen so that the diameter of every Ri, i G {1, . . . , m}'^', is less than 6 and 
Ti-diGjltl] n A) = for all i and j. (Since TCdiA) < oo, there are at most 
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countably many values of t such that 7id{Gj[t] n A) > 0.) Then, we may 
choose 

Vs = {RinA:iG{l,...,mY'} 
and 7 G (0, 1) sufficiently close to 1 such that H40|) holds for Qs = {Qi : i € 
{1, . . . , m}'^ }, where Qi = [jiRi — q) + q) H ^ and q denotes the center of 

For B C A, let 

wb = sup w{x,y) and u;^ = inf w{x,y) 

and define the simple functions ws{x) := YIpgVs^^ ' ^■P(^) ^^'^ ^si^) '■— 
^PePg W-P ■ Xp{x)^ where Xk denotes the characteristic function of a set K. 
Since the distance between any two sets from Qs is strictly positive, Lemma 
Hand equation (|57j) imply 

(41) rsMy^'' > (UQes. q) '"^ > (^Q • 9sAQ)) 

Q^Qg U Q 

Q^H Q&Qs 

Applying Lemma I^J and relation H37() , we similarly have 

(42) fjA)-^i^ < {mp-9jP)y"^' = i-iP-gJP))'"^' 

PeVs A 

Since is a CPD-weight function on ^4 x A, there is some neighborhood G 
of D(A) such that r/ := infgzi; > 0. For 6 > sufficiently small, we have 
P X P C G for all P S Vs, and hence ws{x) > w{x,x) > w^i^x) > rj for 
X £ A. Furthermore, w is continuous at (x, x) G D{A) for "H^-almost all 
X G A and thus, for any such x, it follows that ws{x) and Wsix) converge 
to w{x,x) as 5 ^ 0. Therefore, by the Lebesgue Dominated Convergence 
Theorem, the integrals 

{Ws{x)y^^'' dHdix) and / {wsix))''^^'' dUdix) 



U Q 
QeQs 

both converge to TC^J^{A) as 5 — > 0. Hence, using (|^T|) and we obtain 
(jSHl)- 

Now suppose that Tld{A) > and ujn = {xi , ■ ■ ■ ,xj^} is an asymptoti- 
cally {w, s)-energy minimizing sequence of A^-point configurations on A. It 
is well-known jl81 p. 9] that the weak* convergence result given in H39|) is 
equivalent to the assertion that 
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holds for any almost clopen subset B C A, where we recall that a set B C A 
is called almost clopen (with respect to A and TCd), if the "Hd-measure of the 
relative boundary of B with respect to A equals zero. 

If S C ^ is almost clopen, then the hypotheses of Lemma El imply 

J™ ^^7^ = ^M(^d' (^)) 



Cs,d {nr{B)+nT{A\B))"'''' 



Using relation (|22j) in Lemma [21 and (|38|) for B and A\B, we get 



N^^ N gl^{BYh+g^^^{A\BYh ^ 

showing that ()39() holds. □ 

Theorems 1^ and 121 then follow from Lemma Inland Theorems IT] and IbI as 
we now explain. If s > d and A G is a closed d-rectifiable set, then 
every compact subset B d A\s also closed and d-rectifiable and Theorem^ 
implies that B satisfies condition (|37() and so Theorem [21 then follows from 
Lemmajni If s = d and yl is a compact subset of a d-dimensional C^-manifold 
in , then applying Theorem IbI to every compact subset of A, we get (|S7|) . 
Consequently Theorem follows from Lemma El with s = d. 



5. Proofs of the separation results: Theorem m and 

Corollary [H 

In this subsection we prove Theorem and Corollary ^ For the proof of 
these results we will need Frostman's lemma establishing the existence of a 
non-trivial measure on A satisfying a regularity assumption similar to the 
one in ^Hl for arclength. 

Lemma 7. (see e.g. 'TT, Theorem 8. 8] J. Let a > and A he a Borel set 
in . Then 7{a{A) > if and only if there is a Radon measure /x on M.'^ 
with compact support contained in A such that < fJ-{A) < oo and 

(44) /i[5d/(x,r)] < r", xGR"^', r > 0. 

Moreover, one can find fj, so that fJ,{A) > Cd' .a'^'S' i^) > where Cd'.a > is 
independent of A. 

We proceed using the technique developed in [T6|. Let u;^ := 
{xi, . . . , xat}, G N, > 2, be a (w,s)-energy minimizing configura- 
tion on A (for convenience, we dropped the subscript in writing energy 
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minimizing points Xk,N)- For i = 1, . . . , N let 

TT / \ W{X,Xj) 

Ui{x) := > I 7. , x£A. 

From the minimization property we have that Ui{xi) < Ui{x), x £ A, i = 
1,. . . ,N. If /u is a measure from Lemma[71 set tq := {fi{A)/2N)^^°^ and let 

Di:= A\\jBd'{xj,ro), i = l,...,N. 



Then, by the properties of /i, we have 



/i(A) > KA) - /u [B,, {xj,ro)] > fi{A) - {N - l)r^ > ^ > 0, 



i = 1, . . . , N. Consequently, 
1 



U^{Xi) < 



u^{x)d^,{x) < ^^Y. 



D,, 



w{x, Xj) 



dfi{x) 



Di 



Tdfi{x), i = l,...,N, 



A\Ba,{xj,ro) 



where \\w\\ := sup{|t(;(x, y)| : x,y G A}. Let R := diam A. Then by we 
have fJ-{A) < R^. For every y & A and r E (0, R], using (jl^ we also get 

Ts{y,r) := / -jdfi{x) = / fi{x G A : -j > t}dt 

J \x-y\ J \x-y\ 



A\B^,{y,r) 

^i{A) 



R' 



+ 1^ 



y,t 



-l/s 



s ^a—s 



dt<R''-'+ J t-'^I'dt 



< 



r 



5 



s > a. 



1 + a In s = a. 



N 



Then for i = 1, . . . , N and s > a we have 
(45) 

where Ci > is a constant independent A, w and N. Hence, 

N ,^ II II 



s/a 



=1 
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where Ms^a is a constant independent of A, w, and N. In particular, when 
w = 1, we get 

Since is a CPD-weight function, there are ri,p > such that w{x,y) > rj 
whenever |x — y| < p. Assume that < p and let is and js be such 

that S{uj'^) = \xi^ — Xj^\. Then with some constant C2 > independent of 
N and the choice of we obtain from 1)451) 



> 



Hence, 

where Co = Co (A, ri;, a, s) > 0. Thus, in any case, 

6{uj*n) > mm{p, CoiY^^/"} > CiV"^/", > 2, 

for a sufficiently small constant C^ > independent of N and w^. In 
particular, when w = 1, we have 

The case s = a is handled analogously, which completes the proofs of The- 
orem 0] and Corollary n 

6. Proof of Theorem E] 

The essential ingredient in the proof of Theorem [S] is the following lemma 
which assumes lower regularity. We say that a set K CM."^ is lower a-regular 
if there are positive constants Co and tq so that 

(46) (Co)-V" < w„(i^nSd.(^>0) 

for all X G AT and r < tq. 

Lemma 8. Suppose K C is compact and lower a-regular and a K. 
Further suppose s > a and w : K x K ^ [0, 00] is a CPD-weight function 
on K' X K' for any compact K' C K \ {a}. If w has a zero of order at most 
t at {a, a), where <t < s, then 

(47) glJK) > CiCo-^/"2-(^+*) ^^"(-)) ^^^^ ' 

Proof. Let lon = {xi, ■ ■ ■ , xjy} he a configuration of N distinct points in 
K. For i = 1, . . . , A^, let Pi = \xi — a\, ri = min^yj \xi — Xj|, and choose 
Vi ^ ojn such that \xi — yi\ = ri. Since K is bounded, there is some finite L 
(independent of A^) such that there are at most L — 1 of the points Xi £ ujm 
with the property that > r^. We order the points in wjy so that PN ^ Pi 
for i = 1, . . . ,N and so that < ro for i = 1, . . . , A^ — L. It follows from 
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Cauchy's and Jensen's inequality (see (29) of ^2]) that if 7i,...,7jvf are 
positive numbers, then 

-s/a 

(48) " _ , ^ 



M / M \ 

1=1 \j=i / 



from which we obtain 

N-L , . N-L ^ 

(49) E^M > E > E ^ 

i=i « i=i * 

/N-L ^ \-^/" 

For i = 1, . . . , — 1, observe that 

= min \xi — Xj\ < \xi — a| + min \a — Xj\ < Pi + Pn 1^ ^Pi 

and so 

(50) \x - a\ < \x - Xi\ + \xi - a\ < rj/2 + pi < 2pi, {x G B{xi, rj/2)). 
Using (|46jl and (|5U|) we have 

"""^ < C7o2("/^)("+*) / l^dHaix) 



pf' ~ JKnBi.,,n/2) \x - a|W« 

for i = 1, . . . , N — L. Since B{xi,ri/2) and B{xj,rj/2) are disjoint for i ^ j, 
it follows that 

N-L a f 1 

which combined with (|49p completes the proof. □ 

Remark: If ii' is a-regular at a in the above lemma, then the integral 
|x-a|^('")/'' d7ia{x) appearing in (|^7)) is finite (cf. ^1 p. 109]) and thus 
the Lebesgue Dominated Convergence Theorem gives 

lim / , ,^ , , dTCaix) = 

S-*oJKnB,,{a,5) |x-a| 

and so liuis^o gf^a{K Ci Bii'{a,6)) = oo. 

Now we are prepared to complete the proof of Theorem El First note 
that the hypotheses of Theorem O (namely that A is a^-regular at and 
w has a zero of order of at most t < s at Oj for i = 1, . . . ,n) imply that 
f^w{x,x)-'/'^dnd{x) < oo. 

Suppose e > . By Lemma |H1 and Lemma |21 we can find 6 > such that 
Be ■■= [j7=ii^^ Bd'iai,6)) satisfies Q^J^Be) > (note that if a < d and 
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gl^{K) > 0, then gl^{K) = oo) and ^^(^0 = Ia. wi^.x)-'^'" dUdix) > 
(1- €)n'j'"{A), where A, := A\ B,. 

Since w \s a, CPD-weight function on A^ x it fohows from Theorem [21 
that g'^^{A^) exists and equals Cs^dH'd^ {Af)"^/'^ . Lemmal^then gives 

(51) gl^{A) > {gl,{A,)-<'l^ + gJ^^{B,)-'^l^)-^l'' 

> {c:;'J'nT{A,) + e''/T'''' 

> {C;^J''H'f{A)+e'^/')-'/'^. 

Also, we clearly have 

(52) giM) < alMe) = Cs^nTiA.)-'" < Cm(i - e)-^i'nr{A)-'". 

Taking e — > in and shows that g''^^{A) exists and equals 
Cs,dn'f{A)~'/'^. If Hj'"(A) > 0, then, as in the proof of Theorem 12 
Lemma 121 implies that H15() holds for any asymptotically (w, s)-energy min- 
imizing sequence of configurations cDtv = {^i, • • • , x^}, = 2, 3, . . ., for ^ 
which completes the proof of Theorem El 
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